Quasiparticle relaxation in optically excited high-Q superconducting resonators 
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The quasiparticle relaxation time in superconducting films has been measured as a function of 
temperature using the response of the complex conductivity to photon flux. For tantalum and 
aluminium, chosen for their difference in electron-phonon coupling strength, we find that at high 
temperatures the relaxation time increases with decreasing temperature, as expected for electron- 
phonon interaction. At low temperatures we find in both superconducting materials a saturation 
of the relaxation time, suggesting the presence of a second relaxation channel not due to electron- 
phonon interaction. 
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The equilibrium state of a superconductor at finite 
temperatures consists of the Cooper pair condensate and 
thermally excited quasiparticles. The quasiparticle den- 
sity riqp decreases exponentially with decreasing temper- 
ature. These charge carriers control the high frequency 
(to) response of the superconductor through the complex 
conductivity a± — iu2- At nonzero frequencies the real 
part a i denotes the conductivity by quasiparticles and 
the imaginary part 02 is due to the superconducting con- 
densate [H [3] • When the superconductor is driven out of 
equilibrium it relaxes back to the equilibrium state by the 
redistribution of quasiparticles over energy and by recom- 
bination of quasiparticles to Cooper pairs. The recombi- 
nation is a binary reaction, quasiparticles with opposite 
wavevector and spin combine, and the remaining energy 
is transferred to another excitation. The latter process 
is usually controlled by the material dependent electron- 
phonon interaction 0, 0| . With decreasing temperatures 
the recombination time increases exponentially reflecting 
the reduced availability of quasiparticles. Here, we report 
relaxation time measurements in superconducting films 
far below the critical temperature T c . We find strong 
deviations from exponentially rising behavior, which we 
attribute to the emergence of an additional relaxation 
channel in the superconducting films. 

We have measured the time dependence of the complex 
conductivity of superconducting films after applying an 
optical photon pulse. In addition, the noise spectrum is 
measured in the presence of a continuous photon flux [B[ . 
The superconducting film is patterned as a planar mi- 
crowave resonator. The resonator is formed by a mean- 
dering coplanar waveguide (CPW), with the central line 
3 /im and the slits 2 /Ltm wide, and is coupled to a feed- 
line, see Fig. [T^ 0]. The complex conductivity results 
in a kinetic inductance Lk oc l/duxj?, for thin films with 
thickness d, which is due to the inertia of the Cooper pair 
condensate. It sets together with the length of the central 
line the resonance frequency: luq = 2i:/Al^(L g + Lk)C, 



with I the length of a quarterwave resonator, L g the 
geometric inductance and C the capacitance, both per 
unit length. The variation in kinetic inductance due to 
photons is connected to the quasiparticle density n qp by 
SL)./Lk = \8n qp /n cp , with n cp the Cooper pair density 
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FIG. 1: (a) A quarter wavelength resonator, capacitively cou- 
pled to a feedline, formed by the superconducting film (gray) 
interrupted by slits (black), (b) The resonator exhibits a dip 
in the magnitude and circle in the complex plane (inset) of 
the feedline transmission Six. (c) The feedline transmission 
is converted into a phase 6 and amplitude A using the equi- 
librium resonance circle as reference (right inset). The re- 
sponse to an optical pulse of length 0.5 fjs (at t=0) (open cir- 
cles) exhibits an initial rise due to the response time (3.7 fj,s) 
of the resonator and subsequently follows an exponential de- 
cay (34 (Us) (dashed), reflecting the restoration of equilibrium 
(Eq.[T|. The response is measured with a signal generator, low 
noise amplifier (LNA) and quadrature mixer (upper inset). 
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(n qp <C n cp ). Resonance frequencies used lie between 3-6 
GHz. For a quarterwave resonator at 6 GHz, the length 
of the meandering superconducting CPW-line is 5 mm. 
The resonator is capacitively coupled by placing a part 
parallel to the feedlinc. 

The resonators are made from superconducting mate- 
rials with different electron-phonon interaction strengths, 
tantalum (strong interaction) and aluminium (weak in- 
teraction). The tantalum film, 150 nm thick, is sputtered 
on a high resistivity silicon substrate. A 6 nm thick nio- 
bium seed layer is used to promote the growth of the de- 
sired tantalum alpha phase • The critical temperature 
T c is 4.43 K, the low temperature resistivity p is 8.4 fiilcm 
and the residual resistance ratio (RRR) is 3.0. A 100 nm 
thick aluminium film is sputtered on silicon (T c =1.25, 
p=1.3 fiftcm, RRR=3.7). Alternatively, a film of 250 
nm thick is sputtered on silicon (T c =1.22, p=1.0 pQcm, 
RRR=6.9) and another one of 250 nm is sputtered on A- 
plane sapphire (T c =1.20, p=0.25 pflcm, RRR=ll). The 
samples are patterned using optical lithography, followed 
by wet etching for aluminium and reactive ion etching 
for tantalum. For both materials quality factors in the 
order of 10 6 are reached. The sample is cooled in a cryo- 
stat with an adiabatic demagnetization refrigerator. The 
sample space is surrounded by a cryoperm and a super- 
conducting magnetic shield. Alternatively, the sample is 
cooled in a cryostat with a 3 He sorption cooler without 
magnetic shields. A GaAsP LED (1.9 eV) acts as photon 
source, fibre-optically coupled to the sample box. 

The complex transmission S21 of the circuit is mea- 
sured by sweeping the frequency of the signal applied 
along the feedline (Fig. [1^). Near the resonance fre- 
quency luq the feedline transmission exhibits a decrease 
in magnitude and traces a circle in the complex plane 
(full lines in Fig. [TJa). A non-equilibrium state results 
in a resonance frequency shift and broadening of the dip, 
and a reduction and shift of the resonance circle in the 
complex plane (dashed lines in Fig. Q})). The actual sig- 
nals (filled dot and open circle in Fig. [TJd) are obtained by 
sending a continuous wave at the equilibrium resonance 
frequency luq through the feedline, which is amplified and 
mixed with a copy of the original signal in a quadrature 
mixer, whose output gives the real and imaginary part 
of the feedline transmission (upper inset Fig. [TJ:). The 
non-equilibrium response (open circle), compared to the 
equilibrium response (filled dot), is characterized by a 
changed phase 8 and amplitude A, referred to a shifted 
origin in the complex plane (from the equilibrium posi- 
tion x c ). 

The phase 8 with respect to the resonance circle center 
x c is given by 8 = arctan[Im(S , 2i)/(x c — Re(S > 2i))] and is 
related to the change in resonance frequency by: 8 = 
—4(5^-, with Q the resonator loaded quality factor 
A related change in Lk is given by Suio/uio = —^SLk/Lk, 
with a the ratio of the kinetic to the total inductance. 
The phase 8 is therefore a direct measure of the change 



in complex conductivity (given in the dirty limit by): 

Sa 2 



= -2aQ—(f(E),A 
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with f(E) the electronic distribution function character- 
izing the non-equilibrium and A the superconductor en- 
ergy gap. 

The amplitude A depends predominantly on o\ and 
to a smaller degree on 02- The amplitude is de- 
termined by the complex transmission S21 by: A — 
\/[Rc(52i) — x c ] 2 + Im(52i) 2 / (1 — x c )- On resonance 
S21 = Qc/(Qc + Qu) with Q u ex (J2/01 the unloaded res- 
onator quality factor and Q c the coupling quality factor, 
leading to 
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By measuring A and 8 in the frequency- and time-domain 
we obtain direct information on the relaxation through 
the complex conductivity of the superconducting films. 

A typical pulse response is shown in Fig. [TJ:. The ini- 
tial rise of the phase 8 is due to the response time of 
the resonator. The relaxation shows up as an exponen- 
tial decay. The right inset of Fig. [TJ; shows the evo- 
lution of the response in the transformed polar plane. 
These data are interpreted as governed by one relaxation 
time. This is justified by performing measurements of 
the noise spectrum and applying the analysis by Wil- 
son et al. 5]. Since the superconducting condensate 
and the quasiparticle excitations form a two-level system 
a Lorentzian spectrum is expected, with the relaxation 
time determining the roll-off frequency. If more dominant 
relaxation processes are present, the noise spectrum is no 
longer a single Lorentzian Q . We have studied the super- 
conducting films under exposure to a continuous photon 
flux. Our films are exposed to an optical white noise sig- 
nal due to photon shot noise, resulting in fluctuations in 
f(E). Where a single time r determines the relaxation 
process the phase or amplitude noise spectrum is 
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with P the absorbed power, Ml the photon energy, and 
re t A denoting the responsivity of the phase or amplitude 
to an optical signal. 

The measured noise power spectra of the amplitude 
and phase of a tantalum sample are shown in Fig. [2l In 
equilibrium the amplitude noise spectrum (dashed blue 
line) is flat over the full range, and the phase noise (solid 
blue line) follows l// a with a w 0.25. The amplitude 
noise is due to the amplifier, remaining unchanged at fre- 
quencies far away from wnwhile the phase noise is domi- 
nated by resonator noise 0, , rolling off at a frequency 
corresponding to the resonator response time (0.5 /xs). 
Under a continuous photon flux we observe excess noise 
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FIG. 2: (Color online) The power spectral density of phase 
(solid line) and amplitude (dashed) in equilibrium (blue) and 
under a continuous photon flux (red) at a bath temperature 
of 310 mK. The cross-power spectral density (solid green) un- 
der a continuous photon flux follows a single pole Lorentzian 
spectrum, S oc [1 + (27r/r) 2 ] -1 , with a characteristic time of 
21.7±0.3 fj,s (solid black). The response time of the resonator 
is 0.5 /is. 



in both amplitude (dashed red line) and phase (solid red 
line) that rolls off to the equilibrium value around 8 kHz. 

The difference in noise levels is equal to the difference 
in responsivity: ta/tq = 0.23 (-13 dB), measured for 
this sample. In addition, we estimate, based on 20 pW 
optical power absorbed by the resonator, a phase noise 
level of —94 dBc/Hz due to photon shot noise, which 
is close to the observed value. Thus we conclude that 
the excess noise is due to variations in f(E) induced by 
the photon flux. In order to eliminate the system and 
resonator noise we calculate the phase-amplitude cross- 
power spectral density (solid green line). We find that its 
spectrum is real, indicating that variations in f(E) ap- 
pear as fluctuations in the amplitude and phase without 
relative time delay, and that the data follow a Lorentzian 
spectrum with a single time. The time measured in the 
pulse response (23.0±0.5 (is) agrees with the one deter- 
mined from the noise spectrum (21.7±0.3 /xs). We have 
checked at several bath temperatures and found, also for 
aluminium samples, only a single time. We conclude that 
the relaxation time is the single dominant time in the re- 
covery of equilibrium. 

The measured relaxation times for temperatures down 
to 50 mK are displayed in Fig. [3] The data shown are 
representative for the relaxation times found in all sam- 
ples of different films. In the high temperature regime 
(T/T c > 0.175) the relaxation times increase for decreas- 
ing bath temperature in a similar manner for both tan- 
talum and aluminium samples until a new regime is en- 



tered around T/T c ~ 0.15. The tantalum samples clearly 
show a non-monotonic temperature dependence, exhibit- 
ing a maximum near T/T c ~ 0.15. Two aluminium films 
show a less pronounced non-monotonic temperature de- 
pendence. We do not see a non-monotonic tempera- 
ture dependence in samples of aluminium with the lowest 
level of disorder (highest RRR). Below T/T c ~ 0.1 the 
relaxation times become temperature independent at a 
plateau value of 25-35 /is for Ta, 390 fis for 100 nm thick 
Al on Si, 600 its for 250 nm thick Al on Si and 860 fis for 
250 nm thick Al on sapphire. 

The relaxation times for aluminium are measured in 
half wavelength resonators where the central line is iso- 
lated from the ground plane. For the directly connected 
quarter wavelength resonators a length dependence was 
found. For tantalum the values are found to be length in- 
dependent in both cases. Consequently, the data shown 
are not influenced by quasiparticle outdiffusion. Also, 
the relaxation times remain unchanged when instead of 
an optical pulse a microwave pulse at frequency loq is 
used. In this method only quasiparticle excitations near 
the gap energy are created by the pair-breaking current. 
This observation leads us to believe that the observed 
decay is due to recombination of quasiparticles with en- 
ergies near the gap. 

The exponential temperature dependence for T/T c > 
0.175 is consistent with the theory of recombination by 
electron-phonon interaction jij. The dotted lines in Fig. 
[3] follow the expression for the recombination time, 
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with To a material-specific electron-phonon scattering 
time. We find for 150 nm Ta on Si t = 42 ± 2 ns 
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FIG. 3: The relaxation times as a function of reduced bath 
temperature for 150 nm Ta on Si (■, •), 100 nm Al on Si (a), 
250 nm Al on Si (V) and 250 nm Al on sapphire (0) samples. 
The inset shows the same data on a linear scale. The dotted 
lines are fits to the data using Eq. |4] 
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and for 250 nm Al on Si To = 687 ± 6 ns. The devia- 
tion from the exponential rise and the low temperature 
behavior is incompatible with the established theory for 
electron-phonon relaxation. We assume that an addi- 
tional relaxation channel 11| is dominant at low temper- 
atures, where the electron-phonon mechanism becomes 
too slow. 

In previous experiments using superconducting tunnel 
junctions a similar saturation in the quasiparticle loss has 
been reported. For photon detectors inverse loss rates in 
the order of tens of microseconds have been found for 
tantalum 12, 13, 14, 15 1 and hundreds of microseconds 
for aluminium [6J. Some of these experiments also in- 
dicated a non-monotonic temperature dependence [l6j |. 
Most of these observations have been attributed to trap- 
ping states at surfaces or in dielectrics. The fact that our 
similar experimental results occur in simple supercon- 
ducting films and two different materials suggests that 
processes in the superconducting film itself lead to the 
observed low temperature behavior. 

The observed saturation in the relaxation times in our 
samples is reminiscent of experiments in normal metals 
on inelastic scattering in non-thermal distributions and 
on dephasing in weak localization studies. The apparent 
saturation of the dephasing time and the strong quasi- 
particle energy exchange at low temperatures have been 
shown to be caused by dilute concentrations of magnetic 
impurities [SEE! Sot It is known that in supercon- 
ductors a large density of magnetic impurities decreases 
the critical temperature. For dilute magnetic impurities 
the local properties are most important. In experiments 
with magnetic adatoms impurity bound excitations arise 
2l|, tails in the density of states within the gap might 
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form and the formation of an intragap band with growing 



impurity concentration are predicted [22L |23j . In ongo- 
ing experiments we observe a gradual decrease of the re- 
laxation time with an increasing ion-implanted magnetic 
impurity concentration (0-100 ppm). However, disorder 
plays a role as well and further experiments are needed 
to clarify possible relaxation processes 24|. 

In conclusion, we find that the quasiparticle relaxation 
times, probed by means of the complex conductivity, sat- 
urate for both tantalum and aluminium, below a tenth of 
the critical temperature. We suggest that the saturation 
of the relaxation time is due to the presence of a relax- 
ation channel, which is not caused by the conventional 
process dominated by electron-phonon interaction. 
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